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I. INTRODUCTION 

In 1931 Koopman 0| showed that one can associate with Hamiltonian systems of ordinary 
differential equations the Schrodinger equation in Hilbert space of square integrable functions. 
The observations of Koopman have become an important tool in the spectral theory of dynamical 
systems In 1932 Carleman following ideas of Poincarc and Fredholm demonstrated that 
nonlinear systems of ordinary differential equations with polynomial nonlinearitics can be reduced 
to an infinite system of linear differential equations. This approach is nowadays referred to as 
the Carleman linearization or Carleman embedding. The Carleman approach has been succesfuUy 
applied to the solution of numerous nonlinear problems (see Q and references therein). For 
example it was used for calculating Lyapunov exponents and finding first integrals for the 
Lorenz system In his book Q Varadarajan extended the Koopman linearization to the very 
general case of the phase space replaced with a G-space. When G=R and action of the group 
is given by the flow the Varadarajan observations can be regarded as a generalization of the 
Koopman approach to the case with non-Hamiltonian systems. Recently, such a generalization was 
rediscovered by Alanson [|| who reported the possibility of reformulation of dynamical systems in 
the Hilbert space of square integrable functions. In 1982 Steeb [|| demonstrated that the Carleman 
embedding matrix can be expressed with the help of Bose creation and annihilation operators. 
Inspired by this observation the author introduced in 1987 the Hilbert space approach to nonlinear 
dynamical systems pO| ] . The formalism is based on the reduction of nonlinear dynamical systems 
with analytic nonlinearitics to the abstract, linear Schrodinger-like equation in Hilbert space with 
non-Hermitian boson Hamiltonian. The treatment amounts a far-reaching generalization of the 
Carleman linearization technique which corresponds to the particular case with the occupation 
number representation for the Schrodinger-like equation. On the other hand, it works also in the 
case with partial differential equations. The approach has been developed in a series of papers (see 
monograph ||ll[ and references therein) and it has been shown therein to be an effective tool in 
the study of both ordinary and partial differential equations. 

As remarked by Carleman Q (see also |l^ ) the simplest polynomial ansatz utilized in his 
linearization scheme can be replaced by an orthogonal polynomial. Indeed, the author showed in 
|ll|] that the Koopman linearization is nothing but a version of the Carleman technique with the 
polynomial linearization ansatz coinciding with a multidimensional generalization of the Hermite 
polynomials. It is worthwhile to note that Carleman who discussed the Koopman formalism 
in 1^ and wrote down the linearization transformation coinciding with that given by Hermite 
polynomials did not recognize such interpretation of the Koopman observations. So in both cases 
of the Carleman and Koopman approach the system is linearized via a polynomial ansatz. We 
remark that the polynomials utilized by Carleman are orthogonal ones only in the complex domain. 
Now, using the Hilbert space formalism introduced by the author one can connect the Carleman 
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polynomial linearization ansatz with the Schrodinger-like equation in Hilbert space. The reason of 
the nonhermicity of the corresponding Hamiltonian is the fact that we deal with complex orthogonal 
polynomials. In summary, the common feature of the Carleman and Koopman approaches is an 
orthogonal polynomial ansatz which allows to reduce the nonlinear dynamical systems to the 
Schrodingcr or Schrodinger-like equation in Hilbert space. The aim of this work is to introduce 
a general method for linearization of nonlinear dynamical systems with analytic nonlinearities 
including the Carleman technique and Koopman approach as a special case. Namely, we show 
that the linearization ansatz enabling reduction of dynamical systems to the Schrodinger equation 
is generated by an arbitrary classical orthogonal polynomial. 

The paper is organized as follows. In section 2 we briefly recall the Carleman technique and its 
Hilbert space generalization. Section 3 is devoted to a short exposition of the Koopman lineariza- 
tion. In section 4 we describe the method for linearization of nonlinear dynamical systems in the 
general case of an arbitrary classical orthogonal polynomial. The explicit relations for the concrete 
polynomials are discussed in section 5. 



II. CARLEMAN LINEARIZATION AND THE HILBERT SPACE APPROACH TO 
NONLINEAR DYNAMICAL SYSTEMS 

We first briefly recall the Carleman technique Consider the real analytic system 

x = F(x), (2.1) 
where F : R'^ R'^ and F is analytic in x. Setting 

k 



i=l 

where x = (xi, . . . ,Xk) satisfies (2.1) and rii E Z+ (the set of nonnegative integers), we arrive at 
the linear differential-difference equation such that 

(/)„ = ^ C„in0m- (2.3) 

Since we can introduce an order in the set Z^, therefore (2.3) is equivalent to the infinite linear 
system of ordinary differential equations. In view of (2.2) the finite system (2.1) is embedded into 
the infinite system implied by (2.3). Indeed, it follows from (2.2) that the solution x to (2.1) is 
linked to the solution 0„ of (2.3) by 

Xi ^ 4)e,, i = 1, . . . , fc, (2.4) 

where = (0, . . . , 0, 1^, 0, . . . , 0) is the unit vector of R*^. Therefore, the Carleman linearization is 
also referred to as the Carleman embedding technique. 

We now outline the Hilbert space approach to nonlinear dynamical systems JTo| . Consider the 
system (complex or real) 

z = F(z), (2.5) 

where F : C*^ C'' and F is analytic in z. Let us introduce the vectors in Hilbert space of the 
form 



|z)' = exp > Jz,n |z), (2.6) 

where z fulfil (2.5) and |z) are normalized coherent states (see appendix B). On differentiating 
both sides of (2.6) with respect to time we arrive at the following linear, abstract Schrodinger-like 
equation in Hilbert space satisfied by the vectors |z)': 
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dt 

where M' is the boson operator such that 



-^|z)' = M'|z)', (2.7) 



M' = ^alF,(a). (2.8) 

1=1 

Here aj, aj, i ^ 1, . . . , k, are the standard Bosc creation and annihilation operators, respectively 
(see appendix A). Evidently, the solutions to the nonlinear system (2.5) are linked to the solutions 
to the linear equation (2.7) by 

a|z)' = z|z)'. (2.9) 

It thus appears that the integration of the nonlinear dynamical system (2.5) can be brought down 
to the solution of the linear Schrodinger-like equation in Hilbert space (2.7). 

We now discuss the connection of the Hilbert space formalism with the Carleman linearization. 
Writing the abstract equation (2.7) in the occupation number representation (see appendix A) we 
obtain 

in- ^ M^^z^, (2.10) 
mez*^ 

where Zn = (n|z)' and — (n|il/'|m). Using (2.6) and (B.6) we get 

^, (2.11) 



n 



where z = (zi, . . . , Zk) obeys (2.5). Thus it turns out that the Carleman embedding technique cor- 
responds to the particular occupation number representation within the Hilbert space approach. 
We note that the polynomials (2.11) form the orthonormal complete set in the Fock-Bargmann 
space specified by the inner product (B.7) (see (B.9)). In other words, (2.11) is simply the nor- 
malized version of (2.2). It should be noted that in the case with monomials like (2.11) with 
fc = 1, the passage to the complex domain is necessary for obtaining the complete orthonormal set. 
Indeed, it can be easily demonstrated that there is no Hilbert space spanned by real monomials 
Pn{x) = CnX^ . However, it might be observed that (2.7) and (2.9) hold in the real domain as well. 



III. THE KOOPMAN LINEARIZATION 



This section is devoted to a brief exposition of the Koopman linearization, more precisely its 
generalization mentioned in the introduction, with the help of the Hilbert space approach described 
in ||ri|. Consider the real analytic system 

x = F(x). (3.1) 

Let |x), where x fulfils (3.1), be the normalized eigenvectors of the position operators (see appendix 
C). Using (C.5) and (C.2) we find that these vectors satisfy the linear Schrodinger-like equation 
in Hilbert space of the form 

||x)-M|x), (3.2) 

where the operator M is given by 

fe 

M = -i^p,F,(q), (3.3) 

i=l 
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where pi, qj, i = 1, . . . , k, are the momentum and position operators, respectively. Clearly, the 
solutions to the nonlinear system (3.1) are related to the solutions of the linear equation (3.2) by 

q|x)=x|x). (3.4) 

We have thus shown that the nonlinear dynamical system (3.1) can be cast into the linear 
Schrodinger-like equation (3.2). Moreover, an easy calculation based on (3.2) and (C.3) shows 
that the vectors defined by 



i r* divF(x)(iT| , 
;2 Jo X) 



(3.5) 



obey the Schrodinger equation of the form 



z||x)-=H|x)--, (3.6) 

where the Hamiltonian is 

fe 

H^^Y.^m{<i)+F,{q)p,). (3.7) 

i=l 

In view of (3.5) the solutions to (3.1) and (3.6) are related by 

q|x)~=x|xr. (3.8) 

Thus it turns out that the nonlinear dynamical system (3.1) can be brought down to the abstract 
Schrodinger equation (3.6). Observe that the Hamiltonian H is nothing but the Hermitian part 
of the operator iM. One may ask if we could analogously symmetrize the operator M' given by 
(2.8). The negative answer follows from observations of Kano who showed that whenever 
Hamiltonian is not linear in the Bose creation operators then the coherent states become unstable, 
that is the relation (2.9) is violated during the time evolution. We finally remark that the 
version of the relations (3.5) and (3.6) was originally introduced by Alanson M. 



IV. LINEARIZATION IN HILBERT SPACE AND CLASSICAL ORTHOGONAL 

POLYNOMIALS 

A. Classical orthogonal polynomials 

As we promised in the introduction we now describe a general method for reduction of nonlinear 
dynamical systems to the Schrodinger-like or Schrodinger equation based on the linearization 
ansatz generated by an arbitrary classical orthogonal polynomial. We begin by recalling the basic 
properties of classical orthogonal polynomials. Let Pn{x) be a normalized classical orthogonal 
polynomial, that is we have 

b 

dx W{x)pn{x)p„^{x) = Snm, (4.1) 

where w{x) is the weight function and n, m E Z+. The orthogonal polynomials satisfy the recur- 
rence formula 

Pn+iix) = {AnX + B.n)pn{x) - C„p„_i(x), n = 0, 1, 2, . . . , (4.2) 
where p-i{x) = and 

An = Bn = AniVn+l - r„), Cn = ■ (4.3) 



4 



The coefficients g„ and r„ in (4.3) are given by 

Pn{x) = qnx" + 9',a;""^ + . . . , 

r„ = ^. (4.4) 

We note that the last formula of (4.3) on C„ holds only for the normalized orthogonal polynomials. 
We now specialize to the case with normalized classical orthogonal polynomials. We then have a 
generalized Rodrigues formula jl^ 

1 f7" 

Pn{x) = -——- — [w{x)X% (4.5) 
Knw[x) ax" 

where Kn are constant and X is a polynomial with coefficients independent of n, and the differ- 
entiation formula [Q 

^dpnixl ^ ^ ^X"x)pn{x)+PnPn-lix), (4.6) 

ax ^ 

where 

a„ = nX'(O) - iX'V„, Ar,(in = -Cn{qiK^ + {n-\)X"\. (4.7) 



B. Reduction to the evolution equation in Hilbert space 

We now come to the discussion of the linearization of nonlinear dynamical systems in Hilbert 
space. As it is well-known the orthogonal polynomials form the complete set in [Q. Based 
on this observation we introduce the Herniitian operator x with the complete set of eigenvectors, 
such that 

x\x) = x\x), (4.8) 

and 

{n\x) ^ w{x)ipn{x), (4.9) 

where n G Z_|_, are the basis vectors of the occupation number representation (see appendix A) 
and Pn (x) is a normalized classical orthogonal polynomial. Evidently, the resolution of the identity 
for the states \x) can be written as 

b 

dx\x){x\^I. (4.10) 

On writing the eigenvalue equation (4.8) in the occupation number representation and using (4.9), 
(4.2), (4.3) and (A. 10) we arrive at the following boson realization of the operator x: 

1 \ B 
X = a = H : — (4.11a) 

AN^l^/N An-iVn An 

^ ^ + rjv -rjv+i, (4.11b) 



An^WN An^WN 
where TV is the number operator. Consider now the differential equation 

x^F{x), (4.12 
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where F is analytic in x. Our aim is to study the dynamics of the time-dependent vectors \x), 
where x satisfies (4.12). These states will be seen to be stable with respect to the time evolution 
given by (4.12). that is (4.8) holds at any time. Expanding the time-dependent vector \x) in the 
basis of the vectors \n), differentiating with respect to time and using (4.9), (4.6), (4.2), (4.5) and 
(A. 10) we find that the vectors |a;) satisfy the Schrodinger-like equation in Hilbert space of the 
form 



where the operator M is 



j^\x)=M\x), (4.13) 



M = 7r^. (4.14) 



Here 



TT = - < a 

2 1 ^jv-i 



,[X"{N - 1) + q,K,] - [X"{N - 1) + q,K,] 



-a 



\/N An-i\/N 

+ [X"{N - 1) + qiKi]rN - [X"N + qiKi]rN+i + - X' + 2iVX'(0)l . (4.15) 

J^o ) 

In formulas (4.14) and (4.15) we have not designated for brevity the dependence of X and X' 
on X. It should also be noted that X" is constant for arbitrary classical orthogonal polynomial. 
Therefore, X" is a c- number in (4.15). Now, it can be verified that the following commutation 
relation holds for arbitrary (normalized) classical orthogonal polynomial 

[x,7r]=X. (4.16) 

Therefore, in view of (4.14) we have 

[x,M]=F{x). (4.17) 

Hence, with the use of (4.13) we find that (4.16) is precisely the condition for stability of the 

time-dependent states i.e. (4.8) where x fulfils (4.12), is valid at any time. In other words, the 
solution to (4.12) and the solution of (4.13) are related by (4.8). It thus appears that the solution 
of the nonlinear differential equation (4.12) can be brought down to the solution of the linear, 
evolution, Schrodinger-like equation in Hilbert space (4.13). We remark that (4.16) is equivalent 
to the following identities satisfied by an arbitrary normalized classical orthogonal polynomial: 

[X"{n - 2) + gii^i]r„_i - 2[X"{n - 1) + <zii^i]r„ + [X"n + giKi]r„+i - 2X'{0) = 0, (4.18) 

-j^[X"in - I) + q,K,] - ^[X"(n + i) + q,K,] = X(r„ - r„+i) = X (-^) , (4.19) 

where X(r„ — r„+i) = X{y)\^^^ _^ We finally note that the formula (4.18) is a direct conse- 
quence of the more general relation such that 

[X"{n - 1) + qiKi]rn - [X"n + giXi]r„+i + 2nX'(0) + = 0. (4.20) 



C. Reduction to the Schrodinger equation 

We now demonstrate that (4.12) can be furthermore reduced to the Schrodinger equation. First 
observe that (4.20) is equivalent to the operator formula 

7r + 7r^ = -X'. (4.21) 



6 



Further, owing to (4.16) the Hermitian operator k defined as 

^=i(7r-7rt), (4.22) 

satisfies the commutation relation 

[x, k] = iX. (4.23) 
Using (4.23), (4.22) and (4.21) we find that the Hermitian operator H such that 



2 V ^ ^ , 
is related to the operator M of the form (4.14) by 

ff = i(M + if), (4.25) 

An easy inspection based on (4.13), (4.8) and (4.25) shows that the vectors defined by 

1 r* dF j_ 

|a;)~:=e2 Jo ^"^la;), (4.26) 
where \x) fulfils (4.13), satisfy the Schrodinger equation of the form 

ij^\xr=H\xr. (4.27) 

By virtue of (4.8) and (4.26) we have 

x\x)~= x\x)~ (4.28) 

where x obey (4.12). Thus it turns out that the solution of the nonlinear equation (4.12) can be 
cast into the solution of the Schrodinger equation (4.27). We remark that in view of (4.25) the 
Hamiltonian H coincides with the symmetrization of the operator iM. 

D. Polynomial linearization ansatz 

We note that (4.9) can be regarded as a linearization ansatz for (4.12). Indeed, on writing (4.13) 
in the occupation number representation we arrive at the linear differential- difference equation 
satisfied by {n\x). In order to show that the actual treatment can be regarded as a generalization of 

the Carleman linearization technique we should reformulate it to deal with polynomial linearization 
ansatz instead of the nonpolynomial one given by (4.9). Consider the vectors 

\xy :=w{x)-i\x), (4.29) 

where |a;) fulfils (4.13). Taking into account (4.13) we find that these vectors satisfy the following 
evolution equation in Hilbert space: 

j^\xy = M'\xy, (4.30) 



where M' is given by 



where the operator n' is 



M' = 7r'4^, (4.31) 
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Tt' = -K--—X (4.32a) 

= \ [a ^ X"{N - 1) - [X"{N - 1) + 2q,K,] ^ at 

- qiKi{rN - rjv+i) - j • (4.32b) 



^0 , 

We remark that in view of (4.32a) 

[x,7r']=X, (4.33) 

which leads to 

[x,M']=F{x). (4.34) 
Clearly, the solution to (4.12) is linked to the solution of (4.30) by 

x\x)' = x\x)'. (4.35) 

We have thus shown that the solution to the nonlinear equation (4.12) can be reduced to the 

solution of the linear Schrodinger-Iike equation (4.30). On writing the abstract equation (4.30) in 
the occupation number representation we arrive at the following equation: 

Xn= ^ M^^Xm, (4.36) 
mez+ 

where Xn = {n\x)' and = {n\M'\m). Eqs. (4.29) and (4.9) taken together yield 

Xn=Pn{x). (4.37) 

Thus it turns out that in the particular occupation number representation the presented formalism 
describes linearization of (4.12) by means of the polynomial ansatz (4.37) given by an arbitrary 
classical orthogonal polynomial. 



E. Linearization of multidimensional nonlinear dynamical systems 

We now generalize the actual treatment to the case with multidimensional nonlinear dynamical 

systems. Having in mind the relation (2.2) which is crucial for the Carlcman linearization, it is 
natural to postulate the following multidimensional generalization of classical orthogonal polyno- 
mials: 

k 

Pn{x) =l[pn,{Xi), (4-38) 

i=l 

where Pmixi) is a (normalized) classical orthogonal polynomial and n G Z^^. Evidently, Pn(x) 
form the orthonormal, complete set in the fc-fold tensor product of L^. We can thus introduce the 
Hermitian operators Xi, i = 1, . . . , k, with complete set of eigenvectors |x), x e R*', satisfying 

[xi,xj]=0, i,j = l,...,k, (4.39) 
x|x)=x|x), (4.40) 

(n|x) = (4-41) 

where |n), n e , span the occupation number representation. Clearly, the resolution of the 
identity for the states |x) is 
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j d''x\x){x\= I. (4.42) 

Furthermore, it can be easily checked (see (A.l) and (A. 2)) that the A;-dimensional generahzation 
of (4.11) is 

Xi = tti- 3-7= + 3-7= a| - ^ (4.43a) 

^ + ^—7=4 + - rjVi+1, i = l, ... ,k. (4.43b) 



An.-i^/NI AN,-iVNi 

Consider now the nonlinear dynamical system 

X = F(x), (4.44) 

where F : R'' — ^ R'^ and F is analytic in x. Proceeding as with (4.12) we find that the time- 
dependent vectors |x), where x fulfils (4.44), obey 

^|x)=M|x). (4.45) 

The boson operator M is given by 
where 



I Li—^-y=[X'\N, - 1) + q,K,] - [X"{N, - 1) + q,K,] ^ a] 

+ \X"{Ni - 1) + giKi]rjv. - {X" N, + gii^i]rjv.+i + - X' {x,) + 2iV,X'(0)| , (4.47) 

I = 1, . . . , A;. 

From (4.43), (4.47) and (4.16) (see also (A.l) and (A.2)) it follows that 

\xi,T^j\ = SijX{xj), i, j = 1, ... , k. (4.48) 

Hence 

[x,M] = F(x). (4.49) 

As with (4.16) we find that (4.48) is the condition for the stability of the time-dependent states 
|x). Therefore, the eigenvalue equation (4.40) is valid at any time. In other words, the nonlinear 
dynamical system (4.44) can be brought down to the solution of the linear evolution equation in 
Hilbert space (4.45). 

It is easy to show using observations of section C that (4.44) can be cast into the Schrodinger 
equation. Namely, we introduce the Hermitian operators 

k = U^i-nl), i = l,...,k, (4.50) 

which in view of (4.48) obey 

[Xr,ks] = iSrsX{Xs), T, s = 1, . . . , k. (4-51) 
Taking into account (4.51), (4.46) and the generalization of (4.21) such that 
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*, + *! = -^||^, i=l,.^.,t, (4.52) 

we find that the Hermitian operator H defined as 

i=l 

is finked to the operator M by 

H = i{M +\6xv¥). (4.54) 

Therefore, the vectors 

|x)~:=e^/>'"^''"|x), (4.55) 
where |x) obey (4.45), satisfy the Schrodinger equation 

i||x)~=iJ|x)~. (4.56) 
Obviously, the following relation holds: 

x|x)~= x|x)~, (4.57) 

where x fulfil (4.44). It thus appears that the nonlinear dynamical system (4.44) can be reduced 
to the abstract Schrodinger equation (4.56). 

We finally discuss the multidimensional generalization of the actual treatment in the case with 
the polynomial linearization ansatz (4.37). The resulting formalism generalizes the approach taken 
up by Carleman. Let us introduce the vectors of the form 



\^)' ={^w{x.)-^^ |x). 



(4.58) 



where |x) obey (4.45). Proceeding as in the case with (4.29) we arrive at the following evolution 
equation in Hilbert space satisfied by the vectors (4.58): 

||x)' = M'|x)', (4.59) 

where M' is 

Here the operators tt^, i = 1, . . . , fc, are given by 

K = ^i-\'^X{x,) (4.61a) 
2 w{Xi) 

= I \a,—^-y=X"{N, - 1) - [X"{N, - 1) + 2q,K,] ^ aj 

- qiKi{rN, - rN,+i) - • (4.61b) 



Notice that 
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[xi, iTj] = SijX{xj), i, j = 1, ... , k, (4.62) 
[x,M']=F(x). (4.63) 



Clearly, 

x|x)' = x|x)'. (4.64) 

Thus it turns out that the nonlinear dynamical system (4.44) can be brought down to the lin- 
ear Schrodinger-like equation in Hilbert space (4.59). Writing (4.59) in the occupation number 
representation we get 

*n= ^ M^^x^, (4.65) 

mGZ^i 

where = (n|x)' and = (n|M'|m). Taking into account (4.58) and (4.41) we find 

a;n=Pn(x). (4.66) 

In view of the form of the relations (2.2), (2.3), (4.66), (4.38) and (4.65) it is plausible to treat 
the actual formalism as a generalization of the Carleman embedding technique to the case with 
the linearization ansatz given by an arbitrary classical orthogonal polynomial. We recall that the 
orthogonal polynomials (2.13) are complex. Nevertheless, as we demonstrate in the next section 
the introduced approach works also in such the case. 



V. EXAMPLES 



A. Complex orthogonal polynomials 

In this section we illustrate the general formalism described above by the examples of the con- 
crete classical orthogonal polynomials. Wc begin with the simplest case of the complex polynomials 
(2.13) for A; = 1. As a matter of fact, these polynomials are not considered as classical ones. Never- 
theless, they satisfy the most important relations characteristic for classical orthogonal polynomials 
such as the recurrence, Rodrigucs and differentiation formulas. No wonder that the introduced 
approach covers the case of polynomials (2.13). To see this consider the complex polynomials 

Pn{z) = (5.1) 



Evidently, for these polynomials (see (4.1), (4.2), (4.4) and (2.14)) 

w{z,z*) = e-''', An = ^l=, B„ = 0, Cn=0, (5.2) 

Qn = -p=, Tn = 0. (5.3) 

We remark that the formula (4.3) on C„ does not hold in the case with complex orthogonal 
polynomials. Furthermore, the Rodrigues formula can be written as 

p^(z) = —— --^K^,^*)X"], (5.4) 

where 

X = l, i^„ = (-1)"\/^. (5.5) 
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It can be easily checked with the use of (5.2), (5.3) and (5.5) that the differentiation formula (4.6) 
is also valid by the polynomials (5.1). Now we can write the complex counterpart of relations (4.8) 
and (4.9) of the form 

z\z) = z\z), (5.6) 
{n\z) =w{z,z*)hn{z). (5.7) 

Proceeding as with (4.8) and (4.9) we find 

z = a, (5.8) 

that is the states \z) satisfying (5.6) are nothing but the standard coherent states (see appendix 
B). Referring back to the earlier remark concerning the coefficient C„ in the recurrence formula 
(4.2) we note that the form of the operator z is implied by (5.2) and 

z = a ^a' — . (5.9) 

An-iVN An^/N An 

This formula is a generalization of (4.11) in the case with complex polynomials when the last 
relation of (4.3) does not take place. 
Consider the differential equation 

i = F{z), (5.10) 

where F is analytic in ^. In opposite to the case of classical orthogonal polynomials the time- 
dependent vectors \z), where z fulfil (5.10) do not satisfy the linear evolution equation in Hilbert 
space. Indeed, applying the algorithm described in section 4 we get 

j^\z) = [ReF{z) + a^F{a)]\z). (5.11) 
On the other hand, the vectors defined as 

\zy = w{z,zT''\z) = e'^^^'\z), (5.12) 
which are the complex counterparts of the vectors (4.29) are easily seen to obey 

j^\zy = M'\zy, (5.13) 

where the boson operator M' is 

M' = a'^F{a). (5.14) 

Notice that since X = 1 (see (5.5)), the operator (5.14) has the structure analogous to (4.31). 
Moreover, identifying n' with we see that (4.33) is valid as well. Furthermore, as an immediate 
consequence of (5.14) and (A.l) we find 

[a,M']=F{a). (5.15) 

As with (4.17) we conclude that the eigenvalue equation 

a\zy = z\zy, (5.16) 

where z satisfies (5.10), holds at any time, that is the nonlinear equation (5.10) reduces to the 
linear Schrodinger-like equation in Hilbert space (5.13). 

The generalization of the linearization algorithm in the case with the nonlinear dynamical system 
(2.5) is straightforward. Evidently, the counterparts of (4.38), (4.39), (4.40), (4.41) and (4.58) can 
be written as 
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k 

Zi 



i=l t=l V"*- 

[ai,a-,] = 0, i,j = l,...,k, (5.17b) 
a|z) = z|z), (5.17c) 

(n|z)= |^]J«;(2i,<)5jp„(z)=exp|^-i^|^i|2jp„(z), (5.17d) 

|z)' = |^n^z;(z„<)-ij |z) =exp (^El^^l') l^)" (5.17e) 

Proceeding analogously as in the case of (4.58) we arrive at eq. (2.7). Furthermore, identifying 
TT^ with a| we find that the formula (4.60) takes place also in the case with complex orthogonal 
polynomials. Using (2.8) and (A.l) we get the following counterpart of the relation (4.63): 

[a,M']=F(a). (5.18) 

Clearly, (5.18) leads to the formula (2.9). Finally, the formula (2.11) is a direct consequence of 
(5.17e). We have thus shown that the particular case with the complex orthogonal polynomials 
(5.17a) corresponds to the Hilbert space approach discussed in section 2. 

B. Hermite polynomials 

We now examine the case of the Hermite polynomials within the actual treatment. Let Hnix) 
be the Hermite polynomials. For these polynomials we have (see (4.1)) 

a = —00 , 5 = 00, w{x) = . (5.19) 
Now let Pn{x) be normalized Hermite polynomials, i.e. 

Pn{x)=n-i2-^^Hr,{x). (5.20) 
Vn! 



The polynomials p„ satisfy the recurrence formula (4.2), where 



^" = \/— T' ^" = 0' '^n = J—T- (5.21) 
Vn + 1 

These formulas can be recovered from (4.3) with the help of the relations 

g„ = 7r-i2t^, r„ = 0. (5.22) 
vn! 

The Rodrigues formula for polynomials p„ is given by (4.5), where 

X = l, ii'„ = 7^^(-l)"2t^/;I!. (5.23) 
Using (5.21) we find that the operator (4.11) is 

X = -^{a + a'^), (5.24) 

that is 

x = q, (5.25) 
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where q is the position operator (see appendix C). Furthermore, taking into account (5.21), (5.22) 
and (5.23) we obtain the following formula on the operator (4.15): 



1 

7^' 



n = — (a^ - a). (5.26) 



A look at (5.26) is enough to conclude that 

TT = -ip, (5.27) 



where p is the momentum operator (see appendix C). Therefore, the operator given by (4.22) 
coincides with the momentum operator, i.e. 

k=p. (5.28) 

We now return to the Schrodinger-like equation (4.45). Eqs. (4.46), (4.43), (4.47), (5.23), (5.25) 
and (5.27) taken together lead to the operator M given by (3.3). Thus, the nonUnear system (4.44) 
reduces to the linear Schrodingor-likc equation (4.45), where M is expressed by (3.3). Moreover, 
by virtue of (4.53) and (5.28) the system can be furthermore brought down to the Schrodinger 
equation (4.56), where the Hamiltonian H is given by (3.7). It thus appears that the particular case 
of the Hermite polynomials within the presented approach refers to the Koopman linearization. 



C. Jacobi polynomials 



Our purpose now is to discuss the case of the Jacobi polynomials within the introduced formalism. 



The Jacobi polynomials Pi°''^\x) are specified by 

a=-l, 6=1, w{x) = il-x)°'{l+xf, 
where a > —1 and [3 > —1. Let p„ be the normalized Jacobi polynomials, so 

' {2n + a + P + l)r{n + l)T{n + a + f3 + ly ' 



Pn{x) 



2«+/3+ir(n + a + l)r(n + /?+!) 
The coefficients in the recurrence formula (4.2) are of the form 

2n + a + (3 + 2 f 



Pt^\x). 



Br, 



I {2n + a + p+l){2n + a + p + 3) 
{n + l)(n + a + l)(n + 13 + l)(n + a + P+l) 

{2n + a + l3+ l)(2n + a + P + 3) 
2{2n + a + P)\l {n + l)(n + a + l)(n + (3+ l)(n + a + (3 + l) ' 

2n + a + /3 + 2 I 



2 



n{n + a )(n + /3)(n + a + /3)(2n + a + /? + 3) 



2n + a + P y (n + l)(n + a + l)(n + /3 + l)(n + a + /3 + l)(2n + a + /3 - 1) 
Accordingly, the coefficients g„ and r„ (see (4.4)) are 

2n + a + /3 + 1 



Qn = 



22n+a+/3+ir(n + l)T{n + a + l)T{n + /? + l)r(n + a + f3+l) 

(a - /3)n 



(5.29) 



(5.30) 



(5.31) 



r(2n + a + /3 + l), 

(5.32) 



" 2n + a + P 

The expressions for the polynomial X and the coefficient Kn in the Rodrigues formula (4.5) are 
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X = l-x^, i^„ = (-l)" 



22n+a+/3+ip(^ + l)r(n + a + l)r(n + (3+1) 



(2n + a + l3 + l)r(n + a + P + 1) 
Now taking into account (5.31) we find that the formula (4.11) takes the form 

2 



(5.33) 



X = a 



/ {N + a){N + P){N + a + p) 

2N + a + (3\j {2N + a + (3 - 1){2N + a + + 1) 

2 



(N + a){N + f3){N + a + P) 



2N + a + f3]l {2N + a + l3- 1){2N + a + /3 + 1) {2N + a + l3){2N + a + /3 + 2) 
Furthermore, making use of (4.15), (5.31), (5.32) and (5.33) we get 

2N + a + /3-2 I {N + a){N + P){N + a + P) 



. (5.34) 



+ 



2N + a + P ]] {2N + a + P-l){2N + a + P + l) 
2N + a + f3 + 2 ' 



{N + a){N + p){N + a + p) 



2N + a + P ]] {2N + a + P-l){2N + a + P + l) 
Finally, by virtue of (4.22) and (5.32) we have 



p^-a^ 



{2N + a + p){2N + a + p + 2) ' 

(5.35) 



k = i 



{N + a){N + P){N + a + p) 



+ 



{N + a){N + P){N + a + P) 



{2N + a + P-l){2N + a + P + l) V {2N + a + P- 1){2N + a + P+l) 



(5.36) 



Consider now the system (4.44). Using (4.46) and (5.33) we find that it can be reduced to the 
linear evolution equation in Hilbert space of the form (4.45), with 



(5.37) 



i=l 



where in view of (4.43) and (4.47) Xi and tTj can be obtained immediately from (5.34) and (5.35) 
by formal replacement of a and N by and Ni, respectively. That is f i = x{a = ai, N = Ni) and 

ifi = 7r(a = ai,N = Ni). 

The system (4.44) can be furthermore brought down to the Schrodinger equation (4.56). On 
taking into account (5.36) we arrive at the following form of the corresponding Hamiltonian (4.53): 



H 



F,(i) .mil 



(5.38) 



where ki = k{a = ai,N = N^). 

We now discuss the Hilbert space counterpart (4.59) of the system (4.44). Using (4.32) we get 



2(7V- 1) 



{N + a){N + P){N + a + P) 



+ 



2N + a + p\^ {2N + a + P-l){2N + a + P+l) 

2{N + a + P + l) I {N + a){N + p){N + a + p) 
2N + a + p V {2N + a + p-l){2N + a + p + l) 



2{a- P)N{N + a + P+1) 
{2N + a + p){2N + a + P + 2) " 

(5.39) 



Therefore, the operator M' given by (4.60) takes the form 
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i=l * 

where tt'^ = it' (a = ai,N = Ni). Wc finally recall that Gegcnbauer polynomials, Lcgendre poly- 
nomials and Chebyshev polynomials are special cases of Jacobi polynomials. We also point out 
that the actual treatment cannot be applied in the case with Chebyshev polynomials of first kind 
Tn{x). Indeed, we have T-i{x) = x ^ 0, which implies violating of the recurrence formula (4.2) 
for n = 0. 



D. Laguerre polynomials 

We finally study the case of (generalized) Laguerre polynomials L'^{x) within the actual ap- 
proach. These polynomials correspond to 

a = 0, & = oo, w{x) = x"e-^, (5.41) 

where a > — 1. Let p„ be the normalized Laguerre polynomials, that is 

r(n + i) 



Pnix) 



T{n + a + l) 

The coefficients in the recurrence formula (4.2) are 



Kix). (5.42) 



■ 1 „_ 2n + a + l r — I '"'^^ 



V(n + l)(n + a + l)' v'('^ + + " + 1) ' \ {n^\){n + a^l)' 

(5.43) 

The constants g„ and r„ in (4.3) are of the form 

\JV(n + l)r(n + a + 1) 
The expressions for the polynomial X and the constant Kn in the Rodrigues formula are 

X = x, Kn = ^/T{n+l)T{n + a + l) . (5.45) 
Now, owing to (4.11) and (5.43) we get 

X = -a\/N + a - \/N -\- a a) + 2N + a + l. (5.46) 
Furthermore, making use of (4.15) one obtains 

n = i(aViVT^ - y]vT^ - 1). (5.47) 

Hence, in view of (4.22) 

k = ^{aVN + a - y/N + a a^). (5.48) 

We are now in a position to write down eqs. (4.46) and (4.53) in the case with Laguerre polynomials. 
We have 

k 

M = Y,^i^^ (5.49) 

1=1 
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and 

where Xi = x{a = ai,N = Ni), ni = Tr{a = ai,N = Ni) and ki = k{a = ai,N 
return to (4.30). Using (4.32) we find 

tt' = -VN + a + N. 

Therefore, the operator (4.60) is 

1=1 

where tt^ = w'{a = ai,N = Ni). 

VI. CONCLUSION 

We have introduced in this work the approach generalizing the classical methods for linearization 
of nonlinear dynamical systems such as the Carleman embedding technique and Koopman formal- 
ism. Moreover, in the light of the observations of section 5 the Hilbert space approach developed 
by the author can be also regarded as a special case of the treatment introduced herein. The role 
played by classical orthogonal polynomials in the presented formalism is remarkable. On the one 
hand, the classical orthogonal polynomials have been shown to be the most natural tool for the 
Hilbert space linearization of nonlinear dynamical systems. On the other hand, the stability of clas- 
sical orthogonal polynomials with respect to the nonlinear time-evolution has been demonstrated 
to be one of the properties which actually determine these polynomials via "canonical" algebraic 
relations like (4.16) and (4.21). As a consequence of the general algorithm the new methods have 
been found in this work for linearization of nonlinear dynamical systems connected with an ar- 
bitrary classical orthogonal polynomial, excluding the case with Hermite polynomials which has 
been shown herein to correspond to the classical Koopman approach. We note that in opposite to 
the existing approaches mentioned above, these methods cover the case of systems with the phase 
space different from the whole R'^. For example, in the case when the Hilbert space coimtcrpart 
of the system (4.44) is (4.45) with M given by (5.37), we have the restrictive condition Xi ^ 1, 
i = 1, . . . , k, where x satisfies (4.44). Analogously, (5.49) leads to the requirement that Xj ^ 0, 
i = 1, . . . , k, where x fulfils (4.44). It is suggested that the phase space of the linearized nonlinear 
dynamical system (4.44) should coincide with a subset of (a, 6)*^ (see (4.1) and (4.38)), where (a, b) 
is the interval associated with the corresponding classical orthogonal polynomial. Therefore, the 
case of (5.37) refers to the system (4.44) such that the solution remains in a fc-dimensional cube 
(—1,1)*'. Clearly, an arbitrary system performing finite motion can be reduced to such one by ap- 
propriate rescaling. On the other hand, (5.49) corresponds to systems (4.44) with x e R^, where 
R-i- is the set of positive real numbers. The well-known example of such systems are rate equations 
of chemical kinetics. The experience with the Koopman linearization, the Carleman embedding 
technique and the Hilbert space approach indicates that the formalism introduced herein would 
be a useful tool in the theory of nonlinear dynamical systems. Last but not least, we point out 
that results of this paper would be of importance in the theory of orthogonal polynomials. Indeed, 
it might be observed that the formulas (4.8) and (4.9) as well as the relations (4.16) and (4.21) 
mentioned above seem to provide a new algebraic approach to the theory of classical orthogonal 
polynomials. We also point out that relations (4.51) amount a generalization of the Heisenberg 
algebra (C.l). On the other hand, we have shown earlier (see section 4) that this algebra can be 
interpreted as a condition for stability of solutions to (3.2). 



(5.50) 
= Ni). We now 
(5.51) 

(5.52) 
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APPENDIX A: BOSE OPERATORS AND OCCUPATION NUMBER 

REPRESENTATION 

We recall the basic properties of the standard Bose operators and the occupation number rep- 
resentation. The Bose creation (a''') and annihilation (a) operators, where a''' = (o|, . . . ,0^.) and 
a = (ai, . . . , afc), satisfy the Heisenberg-Weyl algebra 

[ai, ajj] = 6ij, [at, aj] = [a], a]] =0, i, j = 1, . . . , k. (A.l) 

The Hermitian operators Ni = aja,, i = 1, . . . , k, are called the number operators. These operators 
obey 

[Ni, Nj] = 0, [Ni, aj] = -6ijai, [Ni, a]] = Sija], i,j = l,...,k. (A.2) 

Let us assume that there exists in the Hilbert space of states H, a unique, normalized vector |0) 
(vacuum vector) such that 

a|0) = 0. (A.3) 

We also assume that there is no nontrivial closed subspace of H which is invariant under the action 
of the Bose operators. The state vectors |n), n e Z^, defined as 

I") =(114^1 10)' (A-4) 



are the common eigenvectors of the number operators, i.e. 

N|n) = n|n). (A.5) 

These vectors form the orthonormal basis of H, that is 

fe 



(n|m) = [](5„.„., (A.6) 

i=l 

^|n)(„|=/. (A.7) 



nez^ 



The action of the Bose operators on the vectors |n) has the following form: 

ai\n) = ^/ni\n- Bi), at|n) = ^/m + 1 |n + ej), i = 1, . . . , k. (A.8) 

We finally write down the following formulas corresponding to the case with k — 1, which were 
frequently used throughout this work: 

fiN)a = af{N-l), f{N) = fiN - l)a^ , (A.9) 

where A'^ = a^a is the number operator, 

a|n) = \/n|n- 1), a'^\n) = Vn + 1 |n + 1). (A.IO) 
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APPENDIX B: COHERENT STATES 



We now outline the main facts about the standard coherent states. Consider the coherent states 
|z), where z e C*^, that is the eigenvectors of the Bose annihilation operators 



a|z) = z|z). 

The normalized coherent states can be defined as 



z) = exp ^-^ J2 l^^l^^ ^^'^P (j2 1°)' 



(B.l) 



(B.2) 



where |0) is the vacuum vector. The coherent states are not orthogonal. We have 



(z|w) = exp 



k 



These states form the complete (overcomplete) set. Namely, 

/ dMz)|z)(z| =/, 



where 



k ^ 

d/x(z) = TT —d{ReZi) d{Imzi). 

X X 77- 



(B.3) 



(B.4) 



(B.5) 



The passage from the occupation number representation to the coherent state representation is 
given by the kernel 



(B.6) 



On taking into account (B.4), (A. 7) and (B.6) we arrive at the Fock-Bargmann space of analytic 
(entire) functions specified by the inner product 



{4>\^) = j^Jl,{7.) expU^lz.M (^(z*))*^(z*), 



(B.7) 



where (^(z*) ~ (z|0) exp l-^ip) and z* = {zl, . . . , zl). We remark that in view of (B.4) 

and (B.6) the polynomials 



P. 



= n 



i=l 



/ml 



{B.. 



form the orthonormal complete set in the Fock-Bargmann space that is, we have 

i(z ) = ^nm) 



J ^ dfl{z) exp ^- ^ |z,p j (pn(z*))*Pn 



(B.9) 



where Snm = Yii^l ^rumi- 
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APPENDIX C: POSITION AND MOMENTUM OPERATORS 

We finally collect some basic properties of position and momentum operators. The position (q) 
and momentum (p) operators, where q = {qi, . . . ,qk) and p = (pi, . . . ,pk) satisfy the Heisenberg 

algebra 

[qr,Ps] = iSrs, [qr,qs] = \Pr,Ps]=0, V, S = 1, . . . , k. (C.l) 

These operators are related to the Bose operators by 

q=^(a + at), p=-^(at-a), 

a=^(q + ip), at = -i=(q-ip). (C.2) 
An immediate consequence of (C.l) is 

[p,/(q)] = -*^- (C.3) 

Consider now the eigenvectors |q), q G R*^, of the position operators 

q|q)=q|q). (C.4) 
The normalized eigenvectors can be expressed by 

|q) = TT-i exp(iq2)exp[-i(at - V2<if]\0), (C.5) 
where |0) is the vacuum vector. The states |q) form the orthogonal and complete set, namely 

(q|q') = <5(q-q'), (C.6) 
Jd'q\q){<l\=I. (C.7) 

The passage from the coordinate representation spanned by the vectors |q) to the occupation 
number representation is given by the kernel 

(n|q)= (np„.(g,))exp(-iq2), (^.8) 



where Pn{q) are normalized Hermite polynomials, that is 



Pn{q)=n-i2-^^Hn{q). (C.9) 
vn! 
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